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Abstract 

We study the charmless decays B — > AAh, where h stands for vr"'", K~^, K^,K*~^, or K*^, using 
a 605 fb^^ data sample collected at the T{4S) resonance with the Belle detector at the KEKB 
asymmetric energy e^e~ collider. We observe AAK^ and B^ AAK*^ with branching 

fractions of (4.76l[];^^(stat.) ± 0.61(syst.)) x 10"^ and (2.46l[j;^^ ± 0.34) x 10-^ respectively. The 
significances of these signals in the threshold-mass enhanced mass region are 12.4cr and 9.3a, 
respectively. We also update the branching fraction B{B+ AAK+) = (3.38lo:36 ± 0.41) x 10"^ 
with better accuracy, and report the following measurement or 90% confidence level upper limit 
in the threshold-mass-enhanced region: B{B+ AAK*+) = {2.19tlil ^ 0-33) x with 3.7a 
significance; B{B~^ — > AAvr"^) < 0.94 x 10~^. A related search for B^ AAD^ yields a branching 
fraction B{B° AAD^) = (1.05l[];|J±0.14) x 10"^. This may be compared with the large, ~ lO""^, 
branching fraction observed for B^ — > ppW^. The enhancements near threshold and related 
angular distributions for the observed modes are also reported. 
PACS: 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 

The b s penguin loop process plays an important role in rare B meson decays It 
could be sensitive to new physics beyond the standard model due to additional contributions 
from as yet-unknown heavy virtual particles in the loop. Recently the study of the penguin 
dominated baryonic B decays B^ ppK^ ^ and 5° pA.TX~ js] gave intriguing results. 



The proton polar angular distributions in the baryon-antibaryon helicity frame disagree with 
the expectations for short distance h s weak decays {4]. However, in 5 — * ppK* decays j^, 
the K*^ seems to be fully polarized in the hehcity zero state in agreement with the h ^ s 
weak decay hypothesis. The theoretical hierarchies, B{B^ ppK^) > B{B^ ppK*^) 
and B{B^ ppK*^) > B{B^ ppK*^), from the pole model 6|] are experimentally 
established although the predicted branching fraction B{B^ — > ppK*^) is about a factor of 
20 smaller than the experimental measurement. It is therefore interesting to study the 
corresponding branching fractions for B KKK'^*'> decays, the counterparts with protons 
replaced by A's. 

In this paper, we study the charmless three-body decays B AAh, where h stands for 
7r+, K^, K*^, or K*^ I?]. The mode B+ AAK+ has been previously observed 8] 
and presumably proceeds through a 6 — sss process. This decay process can be related to 
B^ ppK^ as shown in Fig.[T](a) and Fig. (U^b). One can simply replace the ud — ud diquark 
pair with an sd—sd pair to establish a one-to-one correspondence between B — > pph and AAh 
decays. A common feature of these decays is that the baryon-antibaryon mass spectra peak 
near threshold as conjectured in Refs. 0, ^j. The meson carries the energetic s quark 
from the h ^ s transition so that a threshold enhancement of the baryon and antibaryon 
system is naturally formed. However, there is another possibility shown in Fig. (U^c) and 
Fig. [T](d), where the A (instead of the K^) carries the s from the b s transition. It is 
interesting to know the role of this s quark in i? ^ AAK^*'' weak decays. 

Since the branching fractions of B ^ AAK and AAtt^ decays are theoretically expected 
at a level 1^ that is detectable with our present data sample, we attempt to determine the 
branching fractions of the various B AAh decays and compare with the latest measure- 
ments for B pph. We also examine the low mass M^^ enhancements near threshold and 
the related angular distributions in order to investigate the underlying dynamics. 
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FIG. 1: Comparisons of possible decay diagrams between — > ppK~^ /B^ pKir and — > 
AAK+. 

II. EVENT SELECTION AND RECONSTRUCTION 
A. Data Samples and the Belle Detector 

For this study, we use a 605 fb~^ data sample, consisting of 657 x 10^5-8 pairs, collected 
with the Belle detector on the T(4S) resonance at the KEKB asymmetric energy e^e~ (3.5 



and 8 GeV) collider Hj. The Belle detector is a large-solid-angle magnetic spectrometer 
that consists of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and an electromagnetic calorimeter composed of CsI(Tl) 
crystals located inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. 



An iron flux-return located outside of the coil is instrumentec 



identify muons. The detector is described in detail elsewhere 



to detect Kl mesons and to 



3. 
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B. Selection Criteria 



The event selection criteria are based on information obtained from the tracking system 
(SVD and CDC) and the hadron identification system (CDC, ACC, and TOF). All charged 
tracks not associated with long lived particles are required to satisfy track quality criteria 
based on track impact parameters relative to the interaction point (IP). The deviations of 
charged tracks from the IP position are required to be within ±0.3 cm in the transverse 
(x-y) plane, and within ±3 cm in the z direction, where the z axis is defined to be the 
direction opposite to the positron beam. For each track, the likelihood values Lp, L^, 
and for the proton, kaon, or pion hypotheses, respectively, are determined from the 
information provided by the hadron identification system. A track is identified as a kaon 
if Lk/{Lk + Lt,) > 0.6, or as a pion if Lt,/{Lk + L.„) > 0.6. This selection is about 86% 
(93%) efficient for kaons (pions) while removing about 96% (94%) of pions (kaons). 
candidates are reconstructed from pairs of oppositely charged tracks (both treated as pions) 
with an invariant mass in the range 485 MeV/c^ < M^+^- < 510 MeV/c^. The dipion 
candidate must have a displaced vertex and flight direction consistent with a originating 
from the interaction point. We use the selected kaons and pions to form K*^ (— i^gTr"*") 
and K*^ (— * K'^tt") candidates. Events with a K* candidate mass between 0.6 GeV/c^ and 
1.2 GeV/c^ are used for further analysis. Similarly, we select A baryons by applying the 
K'^ vertex displacement and flight direction selection criteria to pairs of oppositely charged 
tracks — treated as a proton and negative pion — whose mass is consistent with the nominal 
A baryon mass, 1.111 GeV/c^ < Mp^- < 1.121 GeV/c^ [l3|. The proton-like daughter is 
required to satisfy Lp/{Lp + L^r) > 0.6. This selection is about 97% (95%) efficient for 
protons (anti-protons) while removing about 99% of pions. 



C. B Meson Reconstruction 

Candidate B mesons are reconstructed in the 5+ AAK+, B+ AA7r+, 5° AAK°, 
AAK*^ and B^ AAK*^ modes. We use two kinematic variables in the center 
of mass (CM) frame to identify the reconstructed B meson candidates: the beam energy 
constrained mass Mbc = \/-^beam ~ P\i ^"^^ energy difference AE = Eb — -Ebcam, where 
-E'beam IS the beam energy, and and Eb are the momentum and energy, respectively, of the 
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reconstructed B meson. The candidate region is defined as 5.2 GeV/c^ < Mbc < 5.3 GeV/c^ 
and —0.1 GeV < AE < 0.3 GeV. The lower bound in AE for candidate events is chosen 
to exclude possible background from baryonic B decays with higher multiplicities. From a 



GEANT IJ] based Monte Carlo (MC) simulation, the signal peaks in a signal box defined 
by the requirements 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and |A^| < 0.05 GeV. To ensure 
that the decay process be genuinely charmless, we apply a charm veto. The regions 2.850 
GeV/c2 < Maa < 3.128 GeV/c^ and 3.315 GeV/c^ < M^a < 3.735 GeV/c^ are excluded to 
remove background from modes with rjc, J/ip and Xco? Xci mesons, respectively. According 
to a study of a rare B decay MC sample, the backgrounds in all candidate regions due 
to self cross-feeds (e.g between B^ AAK~^ and B^ AAi^*"^) or due to other rare 
decays such as B^ — > pAir^ , etc., are negligible. The contribution of the B background 
component with S — > 7A is estimated by fitting the AE distribution. This will be included 
in the systematic uncertainty from fitting by comparing the results with and without this 
background component in the fit. 

D. Background Suppression 

After the above selection requirements, the background in the fit region arises dominantly 
from continuum e~^e~ ^ qq {q = u, d, s, c) processes. We suppress the jet-like continuum 



background relative to the more spherical BB signal using a Fisher discriminant 151] • The 
Fisher discriminant is a method that combines n-dimensional variables into one dimension 
by weighting linearly; the coefficients for each variables are optimized to separate signal 
and background. We optimize the coefficients separately in 7 different missing-mass regions 
based on 17 kinematic variables in the CM frame 16|. The missing- mass is determined 
from the rest of the detected particles (treated as charged pions or photons) in the event 
assuming they are decay products of the other B meson. These missing-mass regions are 
defined as <-0.5, -0.5-0.3, 0.3-1.0, 1.0-2.0, 2.0-3.5, 3.5-6.0, >6.0 (GeV/c^). Probability 
density functions (PDFs) for the Fisher discriminant and the cosine of the angle between 
the B flight direction and the beam direction in the T(4S) rest frame are combined to form 
the signal (background) likelihood Cg [Cb)- The signal PDFs are determined using signal 
MC simulation; the background PDFs are obtained from the sideband data: 5.2 GeV/c^ 
< Mbc < 5.26 GeV/c2 or 0.1 < < 0.3 GeV for the KKK+ , AA7r+, KKK^ and KKK*+ 
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modes; 5.23 GeV/c^ < Mbc < 5.26 GeV/c^ or 0.1 < < 0.2 GeV for the AAK*^ mode. 
We require the hkehhood ratio 71 = Cs/{Cs + Ch) to be greater than 0.5, 0.7, 0.3, 0.5 
and 0.65 for the AAi^+, AA7r+, AAK°, AA/r*+ and AAi^*" modes, respectively. These 
selection criteria are determined by optimization of Ug/ y/rig + ni,, where ra^ and rii, denote 
the expected numbers of signal and background events in the signal box, respectively. We 
use the branching fraction ~ 3 x 10"^ (1 x 10^^ for AAtt"*^) in the estimation of and 

use the number of data sideband events to estimate Uf,. If there are multiple B candidates 
in a single event, we select the one with the best TZ value. The fractions of events that 
have muhiple B candidates are 7.4%, 12.6%, 3.7%, 39.4% and 26.8% for the AAi^+, AA7r+, 
AAA''', AAA'*+ and AAi^*° modes, respectively. The systematic errors due to multiple B 
candidates are described later (Sec. IIVD II) . 

III. EXTRACTION OF SIGNAL 

A. Unbinned Extended Likelihood Fits 

We perform an unbinned extended likelihood fit that maximizes the likelihood function 

e-(^AAh+^.?) JL 

L = n (A^AAh^AAh(Mbc„ Ai?,) + iV,,-P,,-(Mbe„ Ai?,)) . (1) 

i=l 

to estimate the signal yields for the AAK^ , AAtt"*" and h.KK^ modes in the candidate region. 
Here PKhh{Pqq) denotes the signal (background) PDF, is the number of events in the fit, 
and N^Xh ^^id Nqq are fit parameters representing the number of signal and background 
yields, respectively. The AAvr"*" mode can contain a non-negligible cross-feed contribution 
from the hAK^ mode, where the is misidentified as a vr"''. Hence we include a KKK^ 
MC cross-feed shape in the fit for the determination of the AAtt"^ yield. The likelihood 
function is more complicated for the AAJ^"* modes, 

L = n (^AAx*^AAi^* + ^'aax.^aax. + Nq^Pm) ■ (2) 

i=l 

since there are contributions from non-resonant B KKKtt decays and one more variable 
in the fit for the Ktt invariant mass, 0.6 GeV/c^ < Mk-k < 1-2 GeV/c^. 
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B. Probability Density Functions 



We take each PDF to be the product of shapes in Mbc and AE (and Mktt, the re- 
constructed invariant mass of kaon and pion, if applicable), which are assumed to be un- 
correlated. Taking B AAK* for example, for the i-th event, P\Xk* 
PAE{A.Ei) X PxT^^MKm)- For the PDFs related to B decays, we use a Gaussian function to 
represent -PA/bc ^^"^ ^ double Gaussian for Pa_e with parameters determined from MC signal 
simulation. The theoretical p-wave Breit-Wigner resonance function is defined by Eqns. [3l 
m and [5l where A is a normalization factor, F is the width of the peak and rriK* , rrix and m^r 
are the nominal masses of the K*, K and vr {2^, respectively. 

rriK* X F X i^f 



where , = ,/ 1 M.L±jg_Zi!4y _ ,„J. (4) 




We use these functions to parameterize the -Pm^^ distributions for K*^ and and use a 



LASS function obtained from the LASS collaboration [171 to model the nonresonant Pm 



distribution. For the continuum background PDFs, we use a parameterization that was 
first used by the ARGUS collaboration /(Mbc) oc x\/l — x^e~^'^^~^^\ to model the 
Pjv/bc distribution with x given by MbcZ-Ebeam and where ^ is a fit parameter. The Pae 
distribution is modeled by a normalized second order polynomial whose coefficients are fit 
parameters. The continuum background PDF for the K* modes, Pmk-„i is modeled by a 
p-wave function and a threshold function, Pm^vt = ''^ x Pp_wavc + (1 — r) x Pthrcshoid and 

^threshold OC {Mk^ -niK - X elciy<(MK.-niK-nin)+C2X{MK.-mK-m^)^] ^J^g^g and C2 

are fit parameters. 
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FIG. 2: Distributions of (with 5.27 GeV/c^ < Mbc < 5.29 GeV/c^) and Mbc (with |A^| < 0.05 
GeV) for (a) ^ AAK+, (b) AAK^ and (c) B+ AA7r+ modes. The dibaryon mass 

M^j^ is required to be less than 2.85 GeV/c^. The solid curves, dotted curves, and dashed curves 
represent the total fit result, fitted signal and fitted background, respectively. The dot-dashed 
curves in plot (c) show the background contribution from the B'^ — > AAK^ mode. 



10 



(a) 




FIG. 3: Distributions of (with 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and 0.816 < Mk-k < 0.976 
GeV/c^), Mbc (with |A£;| < 0.05 GeV and 0.816 GeV/c^ < Mxn < 0.976 GeV/c^) and Mk-k 
(with |A£;| < 0.05 GeV and 5.27 GeV/c^ < A4c < 5.29 GeV/c^) for (a) B+ AAK*+ and (b) 
AAK*^ modes in the threshold-mass-enhanced region. The sohd curves, dotted curves, and 
dashed curves represent the total fit result, fitted signal and fitted background, respectively. 

IV. PHYSICS RESULTS 
A. Fitting Results 

Figures [2] and [3] show the fit results for fi+ AAK+, AAK°, B+ AA7r+, 

AAK*^ and 5° AAK*° in the M^a region below 2.85 GeV/c^ which we refer 
to as the threshold-mass-enhanced region. The resulting signal yields are given in Table [B 
The significance is defined as a/— 21n(Lo / -^^max) ? where Lq and //max are the likelihood values 
returned by the fit with the signal yield fixed to zero and at its best fit value. These values 
include the systematic uncertainty obtained by varying signal PDF parameters by their 1 a 
errors. 
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TABLE I: Signal yields for each decay mode with M^^ < 2.85 GeV/c^. 



Mode 


KKK+ AAK° AA7r+ AAi^:*+ AAK*^ 


Yield 
Significances (a) 


92.7t\li 45.8tli ^■7^nf2 6.54t|i 31.41,^;^ 
16.4 12.4 2.5 3.7 9.3 



B. Observed Branching Fractions 

1. Branching Fractions 

The differential branching fractions as a function of M^^ for the observed modes are 
shown in Fig. |H Tables [III and IIIII give the yields and the corresponding branching fractions 
for each M^^ bin. The yields are obtained from {AE, Mbc(, Mk-k)) unbinned extended 
maximum likelihood fits for each bin of M^^. We find that a threshold enhancement is also 
present for 5° AAA'*° and B° AAK^ decays. We sum the charmless partial branching 
fractions, where the summation excludes bins in the two charmonium regions, to obtain: 
B{B+ -> AAi^+) = (3.38t[!;^^ ± 0.41) x 10-^ B{B^ AAi^O) = (4.76t°:^^ ± 0.61) x 10-^ 
and B{B° AAK*^) = (2.46+°;?^ ± 0.34) x 10"^ For the AAA'*+ mode, we find B{B+ 
AAK*+) = (2.19^o:88 ± 0.33) x 10"^ with 3.7a significance in the threshold - mass-enhanced 
region using the yield in Table [H The differential branching fractions are obtained by 
correcting the yields for the Mj^j;^ dependent efficiency, which is estimated from signal MC. 
Here, we include the efficiency correction for A polarization reported in Ref. j^, 19| as our 



default MC does not include such an effect. The correction factors are 1.17, 1.23, 1.20, 1.22, 
and 1.16 for B+ AAK+, B+ AAn+, B° AAK^, 5° ^ AAK*^ and B+ AAK*+, 
respectively. These factors are obtained in a model independent way. We first use the 
phase space MC sample to obtain the efficiency function in cos6p, where cos6p is the polar 
angle of proton in the A helicity frame. We then use the cos dp distributions in the data 
sideband and signal regions to find their corresponding average efficiencies. With the signal 
yield information from the fit, the model independent signal efficiency can be estimated. 
Fig. [5] shows the differential branching fractions in bins of cos 6p for B^ AAi^+. This 
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TABLE II: Signal yields and branching fractions B (10 ^) in different My^^ regions for B'^ 
and AAK^ decays. The f symbol indicates a charm veto bin. 



AAK+ 



Mj,-^ (GeV/c2 



B+ AAK+ 
Yield B (10-6^ 



5° ^ AA^o 
Yield B (10-6) 



< 2.4 
2.4 - 2.6 
2.6 - 2.85 
2.85 - 3.128 (t) 
3.128 - 3.315 
3.315 - 3.735 (f) 
3.735 - 4.3 
> 4.3 



24.7l:li 



1 65+0-26 

l.DO_Q 23 



17.5 



+5.3 



5.0 



+3.9 
-3.1 



3.7 



1.' 



+3.8 
-3.0 

j+3.3 



20.8 



0.85 



0.61 



+0.20 
-0.18 

+0.18 
-0.16 



18.3 



+5.1 
-4.4 

+4.8 
-4.2 



1.961H? 



o 7Q+0.37 
'""0.35 

14+011 
"-'•-'-^-0.09 

^■^-■--0.11 



5.9 
39.3 



+3.4 



2.7 

f7.0 
-6.3 



0.61 



1.2 



+2.0 
-1.2 



0.08 



0.05 



+0.09 
-0.07 

+0.08 
-0.06 



2 6+2-^ 
^•"-1.8 



-0.35 
-0.28 



o 70+0.67 
' ^-0.60 



o.ii: 

0.26: 
0.22 



0.3l|:^ 0.03^ 



-0.17 
-0.10 

-0.21 
-0.15 

+0.21 
-0.15 

-0.18 
-0.11 



charmless 



103.4 



+12.9 
-11.2 



3.38 



+0.41 
-0.36 



49.1 



+8.6 
-7.0 



4.76 



+0.84 
-0.68 



TABLE III: Signal yields and branching fractions B (10 ^) in different regions for S' 

AAi^*o decay. The f symbol indicates a charm veto bin. 





50 ^ 


AAi^*o 


M^A (GeV/c2) 


Yield 


B (10-6) 


< 2.4 


19-61^:^ 




2.4-2.6 




o.94l°:r, 


2.6 - 2.85 


4.6+1^ 


0.48l°:l 


2.85 - 3.128 (t) 


43.8tf:^ 


A Qc;+0.83 


3.128 - 3.315 




-0.101°:?° 


3.315 - 3.735 (t) 


q c:+3.6 


0.29t°:i 


> 3.735 


7 9+3.6 
' -^-3.1 


-0 61+0-31 

'-'."-'--0.26 


charmless 


nr. q+9.4 


2 46+0-87 

^.^"-0.72 
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FIG. 4: Differential brandling fractions for (a) B+ AAK+, (b) AAK^ and (c) 

AAK*^ modes as a function of M^^. Note that two bins with 2.850 GeV/c^ < M^j^ < 3.128 
GeV/c^ and 3.315 GeV/c^ < M^y^ < 3.735 GeV/c^ contain charmonium events and are excluded 
from the charmless signal yields. The solid histograms are from phase space MC simulation with 
area normalized to the charmless signal yield. 



distribution is not flat but does agree with the theoretical expectation pi 

To verify the branching fraction measurement procedure, we use J/ip K^*^ events with 
J/ip AA in the region 3.070 GeV/c^ < M^a < 3.125 GeV/c^. Using B{J/ip AA) = 



(1.61±0.15)xl0-3 |20j, we obtain branching fractions of (1.30+°|J)x 10^3^ (0.661°;^^) x 10-^ 
and (2.08l[5:42) x 10"^ for B+ J/^K+, 5° ^ J/^fsfO, 5° ^ Jj^K*^, respectively, which 
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FIG. 5: Differential branching fractions us. cos^pforS^ 
region. 



AAK'^ in the threshold-mass-enhanced 



agree with the world average values 20] within 2a including systematic errors that are 
similar to those for the signal mode discussed below. 



2. Polar Angle Distribution 

Figure [6] shows the angular distribution of the A in the AA rest frame for the threshold- 
mass-enhanced region. The yields are obtained from {AE, Mbc) unbinned extended max- 
imum likelihood fits for each bin of cos^^. The angle 6*;^ is defined as the angle between 
the A direction and the direction in the AA pair rest frame. Here, we make a cos^^^ 
dependent efficiency correction and an average correction for the A helicity dependence as 
discussed above. The distribution shows no significant forward peak, in contrast to the 
prominent peak reported in — s> ppK~^ 2]|, which is a unique signature of the intriguing 
result discussed above. 



3. Helicity Distribution 

We study the K*'^ polarization in AAii'*'^ decay, as the K*'^ meson is found to be almost 
100% polarized with a fraction of (101 ± 13 ± 3)% in the helicity zero state in ppK*^ 
decay |5| . To study the polarization, we use a MC simulation to obtain the efficiency as 
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FIG. 6: Differential branching fractions vs. cos 9^ in the AA pair system for — > AAK~^ in the 
threshold-mass-enhanced region. 



a function of cos 6^ in the threshold-mass-enhanced region, where the angle 6^ is defined as 
the angle between the opposite B direction and the direction in the K*^ rest frame. We 
separate the cos 9k distribution into 4 bins for data. We then use {AE, Mbc, Mxn) unbinned 
extended maximum likelihood fits to obtain signal yields in bins of cos 6k and calculate the 
branching fractions for each bin with the corresponding efficiency. Finally, we use Dqq and 
DIq functions, i.e. | cos^ 9k for a pure helicity zero state and | sin^ 9k for a pure helicity 
one (±1) state, to fit this branching fraction distribution. The fit result is shown in Fig. [3 
We find that the K*^ meson is polarized with (60 ± 22 ± 8)% in the helicity zero state. 



4- Upper Limits and Interpretation 

For modes with signal significance less than 4a, we set the corresponding upper limits on 
the decay branching fractions at the 90% confidence level in the threshold-mass-enhanced 



region. Using the methods described in Refs. [2^, |23|], we obtain B{B^ AAK*^) < 
4.98 X 10~^, where the systematic uncertainty has been taken into account. 

Naively, one would expect that the ratio of B{B~^ —>■ AA7r+) to B{B^ —>■ AAK~^) is similar 
to the one of B{B^ ppn^) to B{B'^ ppK^) j^. The anticipated signal yield for B^ — >■ 
AA7r+ is 26.83lt^3- However, we find no significant signal in the threshold-mass-enhanced 
region for B^ AAtt^ and obtain the upper limit B{B^ AA7r+) < 0.94 x 10"^ at the 
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FIG. 7: Differential branching fractions vs. cos 9k in tlie system for B° AAK*° in the 
threshold-mass-enhanced region. The solid curve is the result of the fit. 



90% confidence level. As a cross-check, we measure B{B^ — > AAK'^) = (3.57lj^;54) x 10~^ 
using the misidentified AA/\+ component in the fit. This value agrees well with our 

B{B^ — > KAK^) measurement. 

Our results may indicate that the contribution of the sd — sd popping diagram to B^ —>■ 
AAtt^ shown in Fig. Mjo), is suppressed relative to the ud — ud popping diagram shown 
in Fig. [9](a). In light of this observation, we move to the b ^ c tree diagram (internal W 
emission) dominated decay B^ AAD'^ {us — us popping), which is shown in Fig. [9]|^d). We 
select the 1.852 GeV/c^ < M^+^- < 1.877 GeV/c^ region for candidates and extract the 
B yield. Figure [8] shows the result of the fit. The signal yield is S.SSl^g as with a significance of 
3Aa. The branching fraction ^ AAD^) is (1.05lJ]iI) x 10"^ < 2.60 x 10"^ at the 90% 
confidence level. This branching fraction contrasts with the large, (1.14±0.09) x 10~^ 20|, 
branching fraction observed for B^ ppD^ {uu — uu popping) shown in Fig. [9](c). It appears 
that the diquark pair popping from the vacuum for us — us is considerably suppressed 
compared with uu — uu. 



C. Comparison with predictions and previous measurements 



Table HVl shows a comparison of the branching fractions for B'^ AAA''^, B^ AAtt 



and B^ AAK^ decays to previous measurements 



and to theoretical predictions 



io|, 
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AE (GeV) M^, (GeV/c^) 

FIG. 8: Distributions of AE (with 5.27 GeV/c^ < Afbc < 5.29 GeV/c^) and A/bc (with |A^| < 0.05 
GeV) for the — > AAD^ mode. The result includes the whole region. The solid curves, 

dotted curves, and dashed curves show the total fit result, fitted signal and fitted background, 
respectively. 
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FIG. 9: Possible diagrams that contribute to pp7r+/AA7r+ and ppD^ / AAD^ . 

which are based on the experimental data for B — » pph and B pAh decays. 

The branching fractions of B^ AAK^ and B^ AAtt^ seem to be consistent with 
the theoretical predictions and the results from previous measurements. 
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TABLE IV: Comparison of branching fractions to previous results and theoretical predictions. 



Mode 


Theoretical prediction [10] Previous measurement [8j Results of our study 


AAK+ 
AA^+ 
AAKO 


(2.8 ± 0.2) X 10-6 (2-9lto:TO ± 0.38) x lO"'' {S.SStoH ± 0.41) x 10"^ 
(1.7 ± 0.7) X 10-"^ < 2.8 X 10~6 (< 0.94 x lO"^)^ 
(2.5 ± 0.3) X 10-6 - (4.76l|] ± 0.61) x 10"^ 



f : This value is obtained in the threshold- mass-enhanced region. 



D. Systematic Study 

Systematic uncertainties are determined using high statistics control data samples. 

1. Reconstruction Efficiency 

• Tracking uncertainty: 

Tracking uncertainty is determined with fully and partially reconstructed D* samples. 
It is about 1.3% per charged track. 

• Particle identification uncertainty: 

For proton identification, we use a A — > p7r~ sample, while for K/ vr identification we 
use a D*^ D^7r~^, K~ti^ sample. Note that the average efficiency difference 

for proton identification between data and MC has been corrected to obtain the final 
branching fraction measurements. The corrections are 7.41%, 7.40%, 7.40%, 7.45% 
and 7.48% for the 5+ ^ AAi^+, B+ AA7r+, AAK^, B+ AAK*+ and 

B^ — s> AAK*^{D^) modes, respectively. The uncertainties associated with the particle 
identification corrections are estimated to be 2% for the proton(anti-proton) from A(A) 
and 0.8% for each kaon/pion identification. 

• A Reconstruction: 

We vary the A selection criteria to estimate their impact on the systematic uncertainty. 
The uncertainties from the A mass cut and requirements on kinematic variables are 
1.9% and 1.5%, respectively. For the reconstruction of A and A, we have an additional 
uncertainty of 4.7% in the efficiency for displaced vertex reconstruction. This is de- 
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termined from the difference between A proper time distributions for data and MC 
simulation. 

• Kg Reconstruction: 

Tlie uncertainty in Kg reconstruction is determined from a large sample of £)~ — >• 
Kg7r~ events. We have an additional uncertainty of 4.9% for Kg reconstruction. 

• 71 selection: 

We study the 71 continuum suppression by varying the 7Z cut value from to 0.9 to 
check for a systematic trend. 

• Multiple Candidates: 

The systematic uncertainty in the best B candidate selection is determined by in- 
cluding multicandidate events satisfying the 7Z cut value when obtaining the signal 
yield and the efficiency for each mode. We then take the difference in the branching 
fractions with and without the best candidate selection as the systematic uncertainty. 

• MC statistical uncertainty: 

The MC statistical uncertainty is less than 2%. 



2. Fitting Uncertainty 
• PDF uncertainty: 

A systematic uncertainty in the fit yield is determined by varying the parameters 
of the signal and background PDFs. The assumption of uncorrelated PDFs for Mbc 
and AE is studied by using 2D smoothed histogram functions for both signal and qq 
MC events. The percentage change in the signal yield is about 0.8%. According to 
our MC simulation study, the rare B decays that will significantly affect our signal 
determination are B^ AAK~^ for B^ AAtt^ mode and B AS'^h for all B 
AAh modes. The latter contributes a 0.5% error and is included in the systematic error 
from fitting. The uncertainty in the fit from the M^n PDF for continuum background 
is determined from the difference between the fit results for the B — * ppK* modes using 
analytical functions (a threshold function and a p-wave function) and using the smooth 
function obtained from the Mktt distribution in sideband data jSj. We quote 1% fitting 
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uncertainties for the Mk-k PDF of continuum background in S ^ KKK* modes. We 
quote a 3.2% fitting uncertainty for the Mk-k PDF of non-resonant B hAK-n, which 
is obtained from the difference in the fit results for the ppK-K mode using an analytical 
function(the LASS function) and using a second order polynomial. The second order 
polynomial is (x - Mx,r(iowcrb) ) x {x-MK^^^^ppevh)), where Mi^^(iowcrb) is 0.63325 GeV/c^ 
and Mi^,r(upperb) is 1.8, 2.4 or 3.0 GeV/c^. The uncertainty in the fit for the Mkk PDF 
of non-resonant ppKn is largest (3.2%) when Mx7r(uppcrb) is 1-8 GeV/c^ j5|- The total 
fitting uncertainties for B+ AAK+, B+ AA7r+, AAA'O, 5° AAK*^ and 

B+ AAK*+ modes are 1.5%, 2.0%, 1.5%, 6.0% and 6.0%, respectively. 

• Fitting bias: 

We use 800 simulated MC event sets to measure the difference between the fit result 
and the expected value. The bias is less than 1% for both 2D and 3D fits. 

3. MC modeling 

• Angular distribution of the proton in the A rest frame: 

As described in IIVB 1\ the efficiency uncertainty due to the polarization of A (A) is 
bypassed by using a model independent method based on data. However, to be conser- 
vative, we quote the percentage difference between efficiencies obtained from the model 
independent method and from the theoretically predicted cos6'p distribution [l9|. This 
modeling uncertainty is about 4.3%. 

• Angular distribution of the A in the AA rest frame: 

We choose the most significant mode, B^ AAK^, with M^^ < 2.85 GeV/c^ to 
obtain its cos % HV B 21 distribution, shown in Fig. El Although it deviates significantly 
from a phase space distribution, the overall efficiency difference from a phase space 
MC sample is small since the efficiency versus cos 6/^ is symmetric and flat. We assume 
that this effect is the same for all other decay modes. Thus, the uncertainties from 
the MC modeling for the angular distribution of cos 6/^ are determined to be 0.9%. 

• Angular distribution of kaon in K*'^ rest frame: 

The uncertainties from the MC modeling of the cos 6k angular distribution in the 
AAK*^ mode about 2.5%. This value is determined from the difference between the 
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efficiency in the threshold-mass-enhanced region obtained from the B yields using 
phase space MC event samples and the efficiency calculated from the efficiency dis- 
tribution function, the theoretical PDFs for the K*^ meson and the ratio of the two 
helicity states obtained by fitting to data. 

4- Total systematic errors 

The systematic uncertainties for each decay channel are summarized in Table |Vl These 
uncertainties are summed in quadrature to determine the total systematic uncertainty for 
each mode. 

TABLE V: Contributions to the systematic uncertainty(in %). 



Source 


AAK+ 


AA^+ 


AAK^ 


AAK*+ 


AAK*o 


AAD^ 


Tracking 


6.8 


6.8 


7.8 


9.2 


7.9 


7.9 


Proton ID 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Charged Kaon(Pion) ID 


0.8 


0.8 




0.8 


1.6 


1.6 


A Reconstruction 


4.7 


4.7 


4.7 


4.7 


4.7 


4.7 


A Selection Cut 


2.4 


2.4 


2.4 


2.4 


2.4 


2.4 


K'^ Reconstruction 






4.9 


4.9 






TZ Selection 


5.3 


5.3 


1.3 


3.5 


3.5 


3.5 


Multiple Candidates 


0.9 


0.9 


2.4 


1.4 


1.4 


1.4 


MC Statistics 


1.0 


1.0 


1.5 


2.0 


2.0 


2.0 


PDF uncertainties 


1.7 


2.2 


1.7 


6.1 


6.1 


3.7 


Fitting bias 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


MC modeling 


4.4 


4.4 


4.4 


4.4 


5.1 


4.4 


Secondary decays [20] 


1.6 


1.6 


1.6 


1.6 


1.6 


2.0 


Number of BB Pairs 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


Total 


12.2 


12.3 


12.9 


15.3 


14.0 


12.9 



22 



V. SUMMARY 



Using 657 xlO^BB events, we observe low mass M^^ enhancements near threshold for 
both the AAK^ and AAi^'*'^ modes, with 12 Aa and 9.3o" significance, respectively. We update 
the branching fraction of —>■ AAK^ mode superseding the previous measurement Q], 
and set upper limits on the modes B^ AAK*'^ and B^ AAtt^ in the threshold-mass- 
enhanced region. No significant signal is found in the related mode B^ AAZ)°. All 
the details are summarized in Table IVll The small value of B{B^ AAtt"*"), the large 
value of B{B^ AAK^), and the absence of a peaking feature in the cos6'a distribution 
for B^ — s> AAK^ indicate that the underlying dynamics oi B ^ AAh are quite different 
from those of -B — pph. These results also imply that the s quark from h ^ s penguin 
diagram does not necessarily hadronize to form a K^] the probability of forming a A is 
not negligible. In addition, because E{B^ AAD^{B^ AA7r+)) is much smaller than 
B{B^ —>■ ppD^{B^ —>■ ppTT^)), it appears that diquark pair popping out from the vacuum 
for us — us {sd — sd) is suppressed compared to uu — uu {ud — ud). 
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TABLE VI: Summary of all B AAh results. 



Charmless branching fractions. 


Mode 


Yield fi(10~6) Significances (a) 


B^ AAK*o 
B+ AAK+ 


49.llf;? 4.76lJ]:^8 ± 0-61 12.5 
25.3+?:^ 2.46lg:f^ ± 0-34 9.0 
103.4l|2-9 3.38^5] ± 0.41 16.5 


Results in the threshold-mass-enhanced region. 


Mode 


Yield B{10^^) Significances (a) 


B+ AA7r+ 
B+ AAK*+ 


7.7Qttf2 < 0.94 at 90% C.L. 2.5 
6.54^|37 2.19tlf8±0.33 3.7 
( < 4.98 at 90% C.L. ) 


Related search. 


Mode 


Yield i3(10"^) Significances (a) 


B° AADO 


5.53t||j5 l-05to,u ± 0-14 3.4 
( < 2.60 at 90% C.L. ) 
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